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Abstract. - Zener tunneling in a binary array of coupled optical waveguides with transverse 
index gradient is shown to produce a sequence of regular or irregular beam splitting and beam 
recombination events superimposed to Bloch oscillations. These periodic or aperiodic Zener- 
Bloch oscillations provide a clear and visualizable signature in an optical system of coherent 
multiband dynamics encountered in solid-state or matter-wave systems. 



Light propagation in waveguide arrays has provided in the past few years a conceptually 
and experimentally relevant model to study the optical analogue [1-8] of Bloch oscillations 
(BO) and related effects typically encountered in solid-state or matter- wave systems [9,10]. 
One of the major advantages offered by waveguide arrays as compared to other optical sys- 
tems [11] is the possibility to easily perform a direct visualization of the wave packet dynamics 
in coordinate space [6] which is complementary to most usual frequency-domain (spectral) 
measurements [9,11]. Recent experiments in waveguide arrays [6,8] have reported the beau- 
tiful observation of optical Zener tunneling (ZT), thus motivating the study of multiband 
dynamics in these systems. In presence of a constant field, which yields a BO motion within a 
single band approximation, ZT usually manifests as a cascading of transitions to higher-order 
bands which results in the damping of the oscillatory BO motion or -in the spectral domain- in 
the broadening of the Wannier-Stark (WS) resonances [9,12]. In Ref. [6] a direct visualization 
of decaying BO via ZT to higher-order bands has been reported in a polymer waveguide array 
with a low refractive index contrast, which is analogous to BO decay observed in Bosc-Einstein 
condensates in shallow lattices [10]. In the experiment of Ref. [6], the BO motion of a light 
wave packet in the first band of the array is periodically damped owing to the appearance 
of ZT which manifests itself as regular outbursts of radiation escaping at the turning points 
of the motion where the wave packet reaches the edge of the Brillouin zone in the recipro- 
cal space. This decay of BO via periodic ZT bursts is irreversible and involves a cascading 
of transitions to higher-order bands. However, as discussed in the context of semiconductor 
superlattices [13, 14], a coherent dynamics instead of an irreversible decay process is expected 
when only two bands, coupled by ZT, are involved. In this case the dynamics of occupation 
probabilities in the two bands shows a complex coherent behavior which is related to the 
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Fig. 1 - (a) Effective refractive index modulation (unit cell) of the binary array used in the numerical 
simulations. Parameter values are: a = 20 fim, w\ = 3.2 /im, u>2 = 3.5 /im, An = 0.003, A = 1.44 /im, 
n e ~ 2.138 (note that wi ^ w%). (b) Band diagram for F — 0. n is the band index. The other 
parameter values are the same as in Fig. 1(a). (c) Behavior of interband coupling coefficient 

spacing of the two interleaved WS ladders [14] and may even lead, under special conditions, 
to suppression of ZT. 

In this Letter BO dynamics in a binary array of coupled optical waveguides is theoretically 
analyzed and shown to provide an ideal optical system for a direct observation of two-band co- 
herent Zener-Bloch dynamics, in which ZT manifests itself as an aperiodic or periodic splitting 
and recombination of propagating light beam superimposed to the BO motion. The binary 
array [15] consists of a sequence of alternating single-mode optical waveguides of different 
widths 2w\ and 2u>2, separated by a/2, showing the same refractive index peak change An 
[Fig. 1(a)]. We describe light propagation in the array by a rather standard two-dimensional 
scalar Schrodinger-like equation [16,17], extended to include a transverse index gradient which 
simulates either transverse thermal heating of waveguides [6] or waveguide curvature [4,5]. By 
writing the electric field as E(x, z, t) = ip(x, z) exp(ikn e z — iut) + c.c, where oj is the angular 
frequency of the field, k = lu/cq its wave number in vacuum, and n e an effective index of the 
structure [17], the evolution of the field envelope ip along the propagation z direction satisfies 
the paraxial wave equation 

9-0 h 2 d 2 ip 

ih— = -- — ■— + [n e - n(x)]ip - Fxip, (1) 
oz 2n e ox 2 - 

where n(x)—n e is the effective modulation of the optical refractive index [17] in the transverse 
x direction [with \n(x) — n e \ <C n e ], h = A/(27r) = 1/fc is the reduced wavelength, and F is 
the index gradient coefficient. For circularly-curved waveguides, F is related to the bending 
radius of curvature R by the simple relation F = n e /R [18]. Since the unit cell of the array 
is a dimcr and supports two modes (namely the symmetric and anti-symmetric combinations 
of the modes guided by the narrow and wide waveguides), the band structure shows two 
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Fig. 2 - (a) Gray-scale plot showing beam intensity propagation in a 8-cm long binary array for 
F — 5.52 m , corresponding to a Bloch period zb = 13.04 mm. The input field is a Gaussian beam 
ip(x,0) — exp[— (x + xo)' 2 /wq] with too = 45 /im and xo = 100 /im. The geometrical parameters of the 
binary array are the same as those of Fig.l. (b) Numerically-computed evolution of band occupation 
probabilities Pi and P2 versus propagation distance, (c) Behavior of Pi and P2 as predicted by the 
two-level equations (3) and (4). 



lowest-order bands [15], separated by a gap E g which vanishes for w± — w-i- If the gap is 
large enough, i.e. if w\ is sufficiently different from w-i and the index change An is not too 
small, LZ is usually negligible at low or moderate values of F and BO occurs with a spatial 
period zb = 2nU/(Fa) = X/(Fa) (see, for instance, [4,6,21]). However, for w\ ~ W2, the 
two bands are separated by a small gap and ZT is no more negligible, although it may be 
still disregarded for higher-order bands. A typical band diagram of the binary array [19] 
for F = and w\ ~ wi is depicted in Fig. 1(b) for parameter values typical of Lithium- 
Niobate waveguide arrays [20]. To study the role of ZT on BO dynamics, we numerically 
simulated propagation of a broad Gaussian beam for different values of the refractive index 
gradient F, i.e. of the Bloch period Zb- For the sake of simplicity, normal incidence has 
been assumed, so that at the input plane of the array the first band (n = 1) is mostly 
excited [16]. For very small values of F, a characteristic BO motion with period zb is observed, 
however as F is increased, i.e. as zg decreases, a sequence of beam splitting and beam 
recombination at planes z — zb/2,3zb/2,5zb/2, ■■■ is observed, which is related to ZT as 
discussed below. The sequence is usually irregular, i.e. aperiodic, however at some special 
values of F it shows a nearly-periodic pattern. Since this ZT sequence is superimposed to a 
BO motion with period zb, the resulting beam motion can be referred to as optical Zener- 
Bloch oscillations (ZBO). As an example, Figs. 2(a) and 3(a) show typical numerical results 
demonstrating either nearly-periodic or aperiodic ZBO. In both cases the array length is 8 
cm. Note that, for a curved waveguide array the values of index gradient F used in the 
simulations of Figs. 2(a) and 3(a) correspond to a radius of curvature R ~ 38.7 cm and 
R ~ 44.5 cm, respectively. To relate the observed dynamics with ZT, we adopt a rather 
general technique developed in solid-state physics [12] and expand the wavepacket ip{x, z) as 
a superposition of Bloch states ft) = u n (x, k) exp(ira) of the array, i.e. we set if)(x, z) = 

J2 n f-n/a dKc n (z, K)(p n (x, k), where n is the band index, — ir/a < n < ft/a, u n (x + a, k) = 
u n (x,K), and the normalization condition dx<p^ l ,(x, K?)(p n (x, k) = (<p n '(x, K?)\(p n (x, k)) = 
5 n ,n'8(K — k') holds. The occupation probability of band n is given by P n (z) = J dn\c n (z, k)\ . 
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For the binary array with the band diagram shown in Fig. 1(b) and for near normal beam 
incidence, the main dynamics comprises only the two bands n — 1,2. Figures 2(b) and 
3(b) show the numerically-computed occupation probabilities for the two lowest-order bands 
versus propagation distance, corresponding to the spatial ZBO motion of Figs. 2(a) and 3(a). 
Note that ZT is characterized by sharp and consistent interband transitions at planes z = 
z.b/2, 3zb/2, 5zb/2, ... connected by plateau regions where ZT is small. In these connecting 
regions, the motion of the two wavepackets constructed by the superposition of Bloch modes 
belonging to the two lowest-order bands are independent each other and shows BO. Beam 
splitting and recombination is thus due to abrupt power transfer between the two wavepackets, 
which in the plateau regions move into opposite directions according to the semiclassical 
motion of a Bloch particle. To get deeper insights into the ZBO dynamics and to explain 
the existence of special values of F leading to a periodic motion, we consider the evolution 
equations for the spectral coefficients c n (z, n), which read [12]: 

ih^- = E n (K)c n (z,K) - iF^- - F y^X n iCi(z, k), (2) 

where X n j[t£) — X; „(k)* = (2Tri/a) Jq dx u* n (dui/ dn) and n, I = 1,2 in the two-band ap- 
proximation. The coupling term X\,2 in Eq.(2) accounts for interband transitions, i.e. ZT. 
For a symmetric index profile [n(x) — n(— x)], the Bloch functions u n (x, k) can be taken such 
that X\ t i — X2,2 — and Xi^k) = iO(k), with O real valued and 0(— k) — 9(k). The 
numerically-computed behavior of Xi t 2 for the binary array of Fig. 1(a) is shown in Fig. 1(c). 
For broad beam excitation at normal incidence, the values of coefficients c\ t 2 at the input 
plane z — are approximately given by Ci(0, k) ~ g(n) and C2(0, k) ~ 0, where the spectrum 
g(n) is narrow around n = . In this case, according to the acceleration theorem a solution 
to Eqs.(2) is simply given by C\^{z,k) = fi^(.Fz/K)g{K — Fz/%), where fi^fa) satisfy the 
two-level equations with periodic coefficients 

EiWfr-iFQMh (3) 

E 2 {K)f 2 + iFQ{K)h (4) 

with /i(0) ~ 1 and /2(C)) ~ 0. The occupation probabilities of the two bands are then given 
by P\(z) — \fi\ 2 and Pi{z) — |/2| 2 - ZT corresponds to non-adiabatic transitions at the points 
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Fig. 4 - Behavior of normalized Floquet exponent difference \fj,2 —fii\/a versus Bloch period zb- The 
dotted curve corresponds to periodic ZBO with period 2zb- 

of avoided crossing between the two bands, i.e at k = it /a, Sir /a, 5tt/cl, where the interband 
coupling is largest and the energy separation E%{k) — E\{k) smallest. This explains the 
characteristic dynamics shown in Figs. 2(b) and 3(b), which is rather well reproduced by a 
numerical analysis of the two-level equations [see Figs. 2(c) and 3(c)]. In order to explain the 
existence of periodic ZBO at special values of index gradient F, let us note that, according 
to Floquet theory the solution to Eqs.(3) and (4) is generally not periodic but characterized 
by the two periods 27r/a and 2tt j\[i2 — Mil; where ifii,2(F ) are the Floquet exponents of the 
system. However, for special values of F such that \fi2(F) — [j,i(F)\/a is a fractional number 
N/M (with N and M irreducible integers), the functions Pi^{z) are periodic with period 
Mzb- For instance, for \fi2 — /^.i | = a/2 the evolution of occupation probabilities is periodic 
with spatial period twice the Bloch period zb- Figure 4 shows the numerically-computed 
behavior of \fj,2 — Mi|/ a versus zb = \/(Fa). Note that, at the intersections with the hori- 
zontal dashed line, ZBO are periodic with period 2zb, which is the case shown in Fig. 2. The 
existence of periodic ZBO is closely related to the spectrum of Eq.(l), which for a two-band 
model is given by two interleaved Wannier ladders whose separation depends on the refrac- 
tive index gradient F [14]. We can express the Wannier ladders for Eq.(l) in terms of the 
Floquet exponents ijix t 2 of the two-level equations (3) and (4) by extending the procedure 
of Ref. [12]. We look for solutions of Eqs.(2) in the form c„(z,k) = c n (n) exp(—iEz jfi) and 
impose the periodic boundary conditions c„(0) = c„(27r/a) to find the spectrum E. After 
some straightforward calculations, this yields for the spectrum the two interleaved Wannier 
ladders E„ — ~/ii(F)F + nFa and E m = —fJ,2(F)F + mFa, where m and n are arbitrary 
integers and i//i,2 are the Floquet exponents of the periodic system Eqs.(3) and (4). The 
condition for periodic ZBO is obtained when the spacing F\fi2 — Mi between the two ladders 
is a fractional multiple of the single ladder spacing Fa. 

In conclusion, an optical realization of coherent BO motion superimposed to a periodic or 
aperiodic ZT sequence, leading to ZBO, has been proposed. This coherent dynamics should 
be directly visualized in a binary array with transverse refractive index gradient as a sequence 
of beam splitting and recombination. 
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